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ABSTRACT 

Approximate solutions of the differential equations of 
continuum isentropic gas dynamics may be obtained for 
flows in the equatorial plane of a binary star atmosphere, 
by analogy between gas dynamics and shallow- water flows 
in a rotating water tank. The local inclination of the tank’s 
bottom is proportional to the resulting local gravitational 
force at the corresponding point in the binary-star atmos- 
phere. The analogy requires a constant polytropic coeffi- 
cient of two. This application of the analogy rests on its 
extension to rotating systems with external force fields. A 
theoretical evaluation of the analogy (for the earth’s atmos- 
phere under hydrostatic conditions with polytropic coeffi- 
cient equal to two) yields the correct TT maximum tT height of 
the atmosphere— where density reaches zero. This study 
is motivated by the possibility that it may be applied to the 
early stages of the evolution of the earth-moon system 
after its hypothetical fission. 
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Goddard Space Flight Center 
and 

Burkhard Schulz-Jander* 

Deutsche Versuchsanstalt fur Luft- und Raumfahrt (DVL) 


INTRODUCTION 

According to observation, a binary star system may possess a common atmosphere which en- 
velops both stars. The motion of a particle in such an atmosphere may be treated by the restricted 
three-body problem in celestial mechanics. Several publications have used this approach, e.g., 

Kopal (1956, 1957) and Mrs. Gould (1957, 1959). According to Prendergast (1960), the mean free 
path in these atmospheres is of the order of magnitude of 1 to 10 km; i.e., this characteristic length 
is so much smaller than either the separation of the stars or the radius of one of them that colli- 
sions between particles in the atmosphere must be taken into account. Therefore, Prendergast (1960) 
studied stationary solutions of the differential equations for continuum isentropic gas flows in the 
equatorial plane of the binary system, by neglecting the pressure gradient and the velocity com- 
ponent normal to the ”Lagrangian surfaces.” Subsequently, Huang (1965) employed a steady-state 
celestial-mechanics approach, which neglects the pressure gradient but takes into account colli- 
sions between particles— using the statistical properties of the Jacobian constants of the colliding 
particles. 

The theory presented in this paper was developed because of a desire to study gaseous motions 
in the earth-moon vicinity during the early stages after a fission of their hypothetical parent body. 

By such motions, mass, energy, and angular momentum would be redistributed and possibly ex- 
ported from this binary system. 

The relatively large mean free path quoted above may be assumed to represent an average 
value for the entire atmosphere enveloping the stars revolving around each other. It therefore 
seems desirable to give a rigorous hydrodynamic treatment of the atmosphere, including the pres- 
sure gradient, at least for the regions of maximum density. In addition, the observed eruption 
phenomena in binary systems, e.g., Huang (1965), render a nonstationary treatment of the atmos- 
phere desirable. If these atmospheric motions are always symmetrical with respect to the equatorial 


*Deutsche Versuchsanstalt fur Luft- und Raumfahrt (DVL) West Germany. 
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plane, the differential equations of motion and energy for three-dimensional isentropic gas flows 
reduce rigorously to their two-dimensional versions in this plane; in the continuity equa- 
tion, however, one term must be neglected in this two-dimensional transient approach. Prendergast 
(1960) points out that the total mass of the atmosphere is concentrated about the equatorial plane; 
therefore, limiting the study to flows in this plane appears to be a relatively unimportant restric- 
tion of generality which, also, affects the publications mentioned above. Even in the equatorial 
plane, the time-dependent nonlinear system of differential equations with three independent vari- 
ables precludes any analytic solutions. Numerical solutions reveal general characteristics only 
after very many individual cases have been computed and evaluated. Because of this situation, the 
analogy between gas dynamics and shallow-water motions is employed here to furnish the theoreti- 
cal background for a water-tank experiment simulating transient flows in the equatorial plane of 
the binary-star atmosphere.* This requires a discussion of shallow-water flows and an extension 
of this analogy to flows in rotating systems subjected to external force fields. An analogy such as 
the one presented here should be illustrated by some typical and sufficiently general results. Un- 
fortunately means are not available to carry out these experiments for the time being. Since this 
situation will not change in the foreseeable future, it was decided to publish the paper now. As a 
restricted test of validity and accuracy of the method presented, the hydrostatic density distribution 
in the earth’s atmosphere is obtained in a special application as discussed further on. 


SHALLOW-WATER FLOWS 


| 2 w 



A container (Figure 1) is filled to a certain height with 
water. At time t w = 0, this water tank is rotating about 
the z w -axis with constant angular velocity . The accel- 
eration of gravity g w is assumed to act in the negative In- 
direction. A cylinder with radius F 0w is mounted at the 
bottom of the container; its axis coincides with the z w -axis. 
Along the wetted surfaces of the cylinder and/or the bottom, 
sources or sinks may add water to the system or drain 
it off according to a law that may vary locally and in time. 
(The total amount of water in the container need not nec- 
essarily be constant.) The shape of the bottom will be 
treated later. 


Figure 1— Water tank. , . . . , , , . 

The subscripts w and g denote water and gas, respec- 
tively; the subscript i stands for w as well as g. Dimen- 
sional properties are marked by an upper bar. The angular velocities I., the lengths F 0i , and the 
densities p 0i will be used as constant reference quantities. Non-dimensional Cartesian coordinates 
x., y. , z., velocities u . , v. , w. , times t. , pressures p. , and densities p. are introduced by relations 


*For general information about this analogy, Courant/Friedrichs (1948) and Wehausen/Laitone (I960) are recommended. 
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such as 


Xi r 0i 


Pi = with 


pi 


Po i r o 2 i 


( 1 ) 


^0 w 


5t - and P 0i = P t ( r o g ) • 


A Froude-number 3 . is introduced as the ratio of the centrifugal acceleration r 0i cj . 2 , at distance 
7 0 . from the axis of rotation, to the gravitational acceleration g. : 


3 i = r 0i W i 2/ Bi = const - - (2) 

The symbol D/Dt w denotes the substantial derivative with respect to time t w . According to Lamb 
(1932, p. 4 and p. 318), the equations of motion of the water flow in a system rotating at constant 
angular velocity read as follows: 


Du 

w 


^P w 

Dt 

w 

2v w " X w = 

dx^ 

Dv w 


<?P W 

Dt + 

W 

2u w ■ y, = 

~ ^ 


(3) 


and 


Here Equation 4 was obtained from 


Dw 

w 

DT 


5 Pw 1 
dz 


(4) 


p 

w Dt 


dp 

dz 


(5) 


using Equations 1. The left-hand sides of Equations 3 through 5 express the mass-fixed accelera- 
tions in a system at rest, provided that the quantities of the rotating system are used. According 
to Goldstein (1960, p. 12), the equation of continuity is independent of the chosen coordinate system, 
since it simply states the constancy of a scalar quantity. It thus reads 


d u dv dvs 

WWW 

<9x + dy + dz 

w J w w 


( 6 ) 
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For o; w = 0, the basic simplifying assumptions of the shallow-water theory are (see Wehausen 
and Laitone, I960, p. 667): 

I. Irrespective of flows in the water, the pressure distribution in a vertical column of water 
is the same as in hydrostatics, and 

n. The horizontal velocity components u w and v w are independent of z w . (Considering the 
no-slip condition at the bottom, this assumption can be satisfied only approximately and 
for short times in every version of the analogy.) 

In addition to assuming the validity of I and II for a> w / 0, the following assumption must be 
made: 

in. The hydrostatic pressure distribution holds true even in the presence of motions relative 
to the rotating frame of reference. 

If the fluid rotates as a solid body, its free surface (subscript Ii) is given by 

z ( r ) - const . + T 2 co 2 /2g , 

1 1 \ w / w w/ & w ’ 


where 


(7) 


If any flows relative to the rotating frame of reference (’’disturbances”) are superimposed to this 
solid body rotation, the free surface is given by 

z - z 0 (x,y,t ] ~ in ( x , y , t ) + z TT (x,y ). 

W 2w\w , ' 7 W w / ' w \ w ’ ’ w / II\w J w/ 

If the vertical accelerations Dw w /Dt w caused by these disturbances are negligible because of a suf- 
ficiently small water depth, then Equation 5 is satisfied by the hydrostatic distribution specified in 
assumption D3, 


p w - ^,8, [ z n ( x w- yw) + 7 ?w( x w' y„- ‘J’ 2 ,] + Poo ■ (8) 

where p QQ is the atmospheric pressure acting upon the free surface. Equation 8 has been derived 
by Lamb (1932, p. 318) for rotating shallow water flows. A nondimensional form of Equation 8 is 
given by 


Poo + 


[( 


2w 



(9) 
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The functions F lw = z w -z lw (x w> y w ) = 0 and f 2w = z w - z 2w (x w , y w , t w ) = 0 represent the 
wetted bottom and the free water surface, respectively (see Figure 1). Therefore, 

DF iw/Dt w = 0 . j = 1, 2 , (10) 

according to Lamb (1932, p. 7), since these surfaces are always composed of the same fluid par- 
ticles, and thus 


d z . 
jv 

d t 


+ U . 

J w 


i Z iw 

<?x it 


+ v. 


az jv 


1 , 2 . 


(ID 


Because of assumption II and Equation 11, integration of the continuity Equation 6 with respect to 
z w in the limits z lw and z 2w gives 

<? ( z 2w" z iw) dr. v 1 dr. v l 

~ dt % + L( Z 2w" 2 lw) U wJ + 3^ [( z 2w " Z lw) V wJ = 0 • (12) 


For cT w = 0, Equations 3 and 12 represent the shallow-water theory, which is a first-order 
approximation, by expanding in powers of a parameter e l , defined as the ratio of the squares of 
the characteristic horizontal and vertical scales of length (Wehausen and Laitone, 1960, p. 466 and 
667). Subsequently, the wavelength \ and the amplitude r) w of the disturbance are employed for 
these scales. According to Laitone (1962), a second-order approximation yields in case of A w - oo 
and £ w = 0 


d P w 
d~z w 


a 



and 



where 


Correspondingly, Laitone (1962) shows for\ w - 0 and £> w = 0 



(13) 


(14) 
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Equations 13 and 14 show for = 0 that the following quantities have to be sufficiently small in 
order to allow an inviscid water flow to be approximated by the shallow-water theory: 

a. the hydrostatic depth H w (x w , y w ) = z u - z lw of the water, 

b. the ratio v w /u w of the disturbance amplitude v w to H w , and 

c. the ratio of H /k . 

w/ w 

It is assumed here that conditions a, b, and c also govern the validity of the shallow-water theory 
in the rotating water tank, provided that both and jgrad z lw (x w , y w )| are sufficiently small. 


THE ANALOGY BETWEEN SHALLOW-WATER FLOWS AND GAS FLOWS 

By assuming Zg = Wg = 0, a two-dimensional transient gas flow (subscript g) will be related 
to the three-dimensional transient shallow-water flow (subscript w). By use of h w (1) as defined 
in Equation 13 and taken at station x w 2 + y w 2 = 1, the following functional relationships now are in- 
troduced between the water flows and gas flows: 


x , 


v 


= 


[ z 2w ( x w y„* ‘J- z i«( x ,’ yw)]/ H w( x ' 


(15) 

(16) 


and 




p g = [ z 2 w ( x w y„> O" z iw K- y»)] yss.M 1 ) = 8 23* 


(17) 


The scales of length and time of the two flows under consideration are connected by the quotients 
7„ /?„ and t It = 5 /Z3 = const., respectively. Since the isentropic flow of a chemically homo- 

geneous gas is considered here, there exists afunctional relation p g = Pg (; g ) between pressure 
and density. The special polytropic relation Pg « p« with « = 2, which follows immediately from 
Equations 16 and 17, can only hold true approximately since the kinetic theory of gases shows that 
1< k < 1.67. 


Substitution of Equations 15, 16, and 17 in Equations 3 and 12 yields, by use of Equation 9, 



dp g 

a.) 

V 15 ^ ” 2 v e“ x *j 


P S <3x g 

g 0 \ 

^Pg 

d ( Z l»/ 3 w) 

\jsr, + 2u * _ y g ) 

■ ay g 

P Z dy & 


(18) 


( 19 ) 
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and 


5p g , d . . d 

5T + 53r K u g ) 


( p t V J 


0 . 


( 20 ) 


The gas flow defined by Equations 15, 16, and 17 is two-dimensional, since the third coordinate z w 
of the corresponding three-dimensional water flow represents the density p g of the gas. Because 


it does not follow from Equations 6 and 20 that p is constant, even though this seems to be sug- 
gested by Equations 1 and 15. By the argument following Equation 5, the left-hand sides of Equa- 
tions 18 and 19 are the mass-fixed accelerations in the rotating system of reference, and Equa- 
tion 20 is the two-dimensional equation of continuity of a compressible medium in this system. 

The last terms in the right-hand sides of Equations 18 and 19 may be interpreted as due to a con- 
servative body force, which can be derived from a potential z lw (x gt y g . Thus Equations 18, 19 
and 20 are the equations of motion and continuity (in the continuum regime) of two-dimensional gas 
flows subjected to conservative body forces in a system rotating at a constant angular velocity. 

The first law of thermodynamics yields for any mass element of the supposedly isentropic gas 
flow the relation c pg dT g = dp g /p g between the specific heat c pg at constant pressure, the absolute 
temperature T g , the pressure p g , and the density of the gas. Therefore, the energy equation in 
the rotating system, because of the argument following Equation 5, reads 


DT 

p c — 
g pg Dt 



if 


const . 


( 21 ) 


Equation 21 is satisfied by two-dimensional, transient calorically perfect gas flows considered here. 
Thus the relation between water and gas flows as given by Equations 15, 16 and 17 is a complete 
analogy within the limitations of the assumptions mentioned earlier. 

The generalization of the analogy developed earlier consists of admitting (a) conservative 
body forces, which determine the bottom shape z w = z lw (x w , y w ) of the water tank, and (b) rotating 
systems. Variable bottom configurations have been used in many applications of the analogy to 
problems without body forces in order to admit arbitrary constant values of k. From 1952 to 1962, 
this issue initiated a discussion in which body forces were not mentioned at all. In agreement with 
the derivation presented earlier, Laitone (1952, 1953, 1961, 1962) and Wehausen and Laitone (1960) 
pointed out the necessity of z lw (x w , y w ) = const, in cases without body forces which were the only 
ones they considered in this context. Loh (1959, 1962) and Bryant (1960, 1962) remarked that 
z iw (x w , y w ) ^ const, is admissible in view of the other approximations involved in the analogy, 
such as constant values of 


7 



The conditions for the validity of the shallow-water theory mentioned earlier also apply to the 
analogy. Experiments by Laitone (1952) confirmed that the analogy yields satisfactory agreement 
between gas dynamics and water flows, for tap water with usual surface tension, as long as the 
water depth is less than 1/4 inch and the wavelength \ w of the disturbances exceeds 3 inches; see 
also Laitone (1953, 1961, 1962). Under consideration of the surface tension of water, Gupta's 
theory (1965) yields about the same limits for the validity of the analogy carried out in nonrotating 
water tanks. 

For sufficiently small values of S3 w , the wave systems caused by the tank's rotation (e.g., Lamb, 
1932) presumably are negligible as compared with the wave systems discussed earlier for £ w = 0. 
Since, in addition, a transient boundary layer growth at the tank bottom is caused by the waves dis- 
turbing the solid-body rotation of the water, the relation t - t" w ^ w following from Equations 1 
and 15 prohibits atmospheric motions to be studied in the water tank beyond sufficiently short in- 
tervals t after the beginning of these waves. The required small values of 7 w demand a suffi- 
ciently small maximum inclination of the tank bottom. This may necessitate the addition of a de- 
tergent to the water because of the surface tension. 

APPLICATION TO AN ASTROPHYSICAL AND A GEOPHYSICAL PROBLEM 

The right-hand sides of Equations 18 and 19 were interpreted as representing body forces 
which result from a potential z lw (x g , y g )/3 w . The gravitational potential of spherically symmetric 
celestial bodies will now be employed to demonstrate the analogy. Figure 2 shows a binary system 

of two celestial bodies rotating about their com- 
mon center of gravity s at the constant angular 
velocity <j g . The masses of the bodies are M a 
and M*; 7 a and denote the distance of any 
point Q(x, y, z) from the centers of the two 
masses, respectively. The gravitational poten- 
tial in the space outside M a and then is given 
by 

rM a ru 

'P = TT + 77 ’ (22) 

where r is the universal gravitational constant. 
In the following, the contribution of atmospheric 
matter to the gravitational potentials will be neglected. The mass ratio ^ = m /3/ /m cx is introduced, 
and the radius of any one of the two celestial bodies is used as a constant reference length 7 0g ; 
i.e., x a - 7 a /7 0g and = 7^/7^ . Since w has the units (cm 2 sec" 2 ) the expression 7 0 2 oJ g 2 = const, 
is employed to obtain the nondimensional potential 



^9 
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where the Froude number 


3 


a 


Y oV 


rM 


a 


is the ratio of centrifugal to gravitational forces at the surface of the celestial body M a . In applica- 
tions of the analogy, 3 a = 3 . In the equatorial plane, the components of the gradient of this poten- 
tial will now be equated to the last terms in Equations 18 and 19 to obtain an expression for the 
shape of the bottom of the water tank: 


grad g - grad m 


(24) 


With z 0 as a constant of integration, it is found that 



The analogy as expressed by Equations 15, 16, and 17 enables one to study in the equatorial 
plane two-dimensional atmospheric motions subjected to a gravitational field ^(x & , y g , z g ) by ex- 
perimental simulation in a water tank with the bottom z w - z lw (x w , y w ) . These atmospheric flows 
are restricted to the equatorial plane z =0, where the complete differential equations (18), (19), 
and (21) take on special two-dimensional forms, provided that the motion is isentropic and sym- 
metrical with respect to the plane z =0. Because of this symmetry, the first derivatives with 
respect to z g of u , v , p , p , and T g vanish identically as functions of x g , y g , and t in the 
plane z g = 0. Since <9w g /dz g £ 0 for Zg = 0, the three-dimensional equation of continuity retains 
this term dw g /dz g , which is absent in Equation 20. If 

[( dW J dZ *)] Zg , 0 

is known as a function of x g , y g , and t g by use of some information outside the analogy, this term 
can be represented by a suitable distribution of sources and sinks in the bottom z w = z lw (x w , y w ) 
of the water tank. The equation of continuity (Equation 20) then possesses a corresponding addi- 
tional source term. 

At present, the authors cannot experimentally test the accuracy of the analogy; a theoretical 
check is therefore desirable. Such a test of accuracy is provided by the application of the analogy 
to the earth's atmosphere under conditions of a known density distribution; i.e., in case of the hy- 
drostatic atmosphere. In this case, the foregoing equations simplify somewhat, since m = 0 and 
rM a /r 0 2 g = g 0 , where g 0 is the gravitational acceleration at the surface of the earth. Then the 
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quotient 3 w /3 a reduces to an expression defining the constant A: 


3 w 


tt 


— — 2 
r n co * 
Ow w 

~r co ^ 
Og g 


A , 


(26) 


i.e., the ratio of the centrifugal forces in the container to those at the earth T s equator. Because of 
Equation 15, r a = r g T r w . Thus, by taking z lw = 0 at r w = 1, Equation 25 yields 

z iw = < 27 ) 

The atmosphere enveloping the earth is assumed to rotate as a solid body, such that in the analogy 
z 2 w = z riw 311(1 V* = 6. From Equation 7 it is thus found 

nr 

o 

W 0 

z 1 1 w " const - + T r w * (28) 


Because of Equations 16 and 17, 


p 

A 


g 

g 


2'3 w ( z i i w z iw) • 


(29) 


The equation of state for a thermally perfect gas, p 0g = p 0g R g T 0g , requires at the surface of the 
earth (subscript 0) 


p 0 g _ R g T o g 1 

p 0g &o r 0g 


Because 3 a - 3 g , and z lw = 0 at r w = 1, Equations 29 and 30 yield 


with a constant B defined by 



2 AB , 


B 


R g > 

g 0 r 0g 


From Equations 27, 28, and 31 finally, 


Z TT “ Z . 

1 1 w 1 w 


_2B + ^(r w z -l)-(l-^-) 


(30) 


(31) 


(32) 


(33) 
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is obtained to express the nondimensional density distribution in the earth’s atmosphere under 
hydrostatic conditions and for « = 2. From left to right, the terms between the brackets in Equa- 
tion 33 represent the contributions of the earth’s surface temperature, the centrifugal force due to 
the assumed rigid-body rotation, and the gravitational force. 

Because of the small ratio of the maximum extension of the atmosphere to the earth’s radius, 
it is reasonable to write r g = r w =i+e, where second and higher powers of e are to be neglected. 
From 1/(1 + e) = l-e+e 2 • * * , the following relation is obtained to replace Equation 33: 

z iiw - z iw = A [2B - ( 1 - 3 g ) e] • (34 ) 

Because of Equation 17, the pressure variation in the earth’s atmosphere according to the analogy 
is 


Z I Iw 

~ Z lw 

2 


r 

1-3 1 

g 

( Z II w 

‘ Z tw)o_ 



1 “ 

2B e 


(35) 


where the effect of the centrifugal force is negligible even at the earth’s equator since 3 
= 3.32 x 10 -3 . Because B = 1.27 x 10' 3 for a surface temperature T 0g = 293 (°K), the analogy 
predicts that the atmospheric pressure p g vanishes at a maximum altitude as given by 


H 


max 


Og 


e 

max 


r o g 


2B 1 

T^IT * 394 r o g 

g 


% 16 km 


It will now be shown that this result for H max agrees with the one of the well-known hydrostatic 
polytropic atmosphere with k = 2, which is defined (a) by neglecting rotation (i.e., 3 = 0) and (b) 

by assuming a constant acceleration of gravity. Under these conditions, Prandtl/Tietjens (1957, 
p. 33) show that 


p g _ / x- i jA "- 1 

Po B V~ K \) 

where, in the notation of this paper, 

h = r g - T 0g , and h 0 = Br 0g . 

Because e = h/7 0g , Equation 26 yields, in the case k = 2, the relationship 



(36) 


(37) 


(38) 
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which is identical to Equation 35 for 3 g = 0. It should be recalled here that the polytropic hydro- 
static atmospheres defined by conditions (a) and (b) mentioned earlier possess a finite maximum 
altitude H max , with the exception of k - 1. 

A BRIEF REVIEW OF RESULTS ON MOTIONS IN THE 
EQUATORIAL PLANE OF A BINARY STAR SYSTEM 

Prominence activity in the sun is a well-known phenomenon. Observation of binary stars sug- 
gests that the prominence activity of the secondary star is stronger on the hemisphere facing the 
primary star. Appropriate initial velocity vectors of these prominences could result in a gaseous 
stream toward the primary or in an escape from the binary system (Sahade, 1960). 

Within the context of the restricted three-body problem of celestial mechanics, several authors 
have studied trajectories of ejected particles in the equatorial plane. For a particle moving in this 
plane, the Jacobi integral shows the existence of curves of zero velocity in a frame rotating with 
the binary system. These closed curves are the "Roche equipotentials" on which the geopotential 
is constant, i.e., the sum of the gravitational potentials of the component stars and the potential of 
the centrifugal force caused by the rotation of the system. These Roche equipotentials define ad- 
missible cross-sections of each component star in the equatorial plane. The equipotential enclos- 
ing the largest area in this plane and only one component star is called the Roche limit. For a 
particle having a given value of the Jacobi integral, the corresponding Roche equipotential repre- 
sents a barrier which the particle, if moving inside its folds, can never penetrate. The trajectories 
of particles ejected from a component star completely filling its Roche limit take place in an area 
of the equatorial plane whose Roche equipotential boundary envelops both component stars. For 
this case, trajectories of particles have been computed by numerical integrations of the equations 
of motion of the restricted three-body problem (Kopal, 1959). The results could be summarized 
as follows: (1) particles fall upon the companion star; (2) particles fall back upon the parent star; 
and (3) particles escape from the system if they possess sufficiently high initial velocities, generally 
of the order of several hundred kilometers per second. 

According to Sahade (1960), velocities of 700 km/sec have been observed in the binary system 
V 444 Cygni. Because of the smaller masses, considerably smaller velocities would be sufficient 
for escape of matter from the earth-moon system after fission of its hypothetical parent body. 

These results of celestial mechanics previously referred to do not account for interactions 
between particles. Therefore, Struve (1957), has pointed out that it would be important to consider 
the problem at hand also from a hydrodynamic al point of view and to allow for the effect of radia- 
tion pressure and magnetic forces. The theory presented earlier in this paper furnishes a back- 
ground for a water tank experiment that attempts to simulate gaseous streams in the equatorial 
plane of a binary star system by accounting for gravitational, inertial, and hydrodynamic pressure 
forces. These results would pertain to the continuum flow regime, which is separated from motions 
caused by applications of celestial mechanics by the rarefied-flow region. The transient water tank 
experiment, therefore, would furnish a necessary supplement to published results from (a) applications 
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of celestial mechanics (Kopal, 1959), and (b) stationary continuum flow theory under rather restric- 
tive assumptions (Prendergast, 1960). 

COMMENTS WITH RESPECT TO THE EXPERIMENTAL SIMULATION OF THE FLOW 
IN THE EQUATORIAL PLANE OF THE BINARY STAR SYSTEM 

In the water tank, the primary and the secondary stars are represented by cylinders mounted 
parallel to the table's z-axis of rotation (Figure 3). Since the primary star, by definition, possesses 
the larger mass, the bottom of the tank reaches 

its smallest z -coordinate at its intersection j w 

with the "primary” cylinder. If the surfaces t 

of the cylinders are sufficiently smooth, the 
motion of the water is affected very little by 
rotations of the cylinders about their own axes 
of symmetry. Therefore, the wetted portions 
of these surfaces may have to be roughened 
artificially if this eigen- rotation of the cyl- 
inders is to influence the water flow appreci- 

ably. Even though the gravitational fields ex- Figure 3_ Sch ematic configuration of water tank sug- 

tend to infinity in the equatorial plane, it is gested to study the atmosphere around the celestial 

sufficient to place a board parallel to the axis bodies shown in Figure 2. 
of rotation of the system on the bottom of this 

tank, provided that the board is located beyond the synchronous satellite orbit of the system. 

A water sink may be placed at the intersection of the board and the bottom. 

Under hydrostatic conditions, the water depth at the circumference of one of these cylinders 
may be prescribed arbitrarily. Because of Equations 32 and 33, this depth is proportional to the 
temperature of the atmosphere at the surface of the corresponding component star. This tempera- 
ture should be sufficiently low to ensure the gaseous state of the atmosphere since Equations 16, 

17, and the equation of state of a thermally perfect gas show that the temperature of the atmosphere 
is proportional to its density in the analogy employed here. The tank's bottom between the cylinders 
does not have to be covered by water continuously in time. As far as mass transfer from the sec- 
ondary to the primary star is concerned, it is sufficient to have intermittent water jets emanate 
from the "secondary” cylinder. The mass flow from the secondary to the primary cylinder can be 
visualized by adding dye to the water pouring out of the secondary cylinder. 

CONCLUDING REMARKS 

The objective of conducting a water tank experiment should be a qualitative simulation of 
transient astrophysical effects such as those mentioned earlier in this paper. The following de- 
tails militate against quantitative validity of the results: 

a. the assumption of a constant polytropic coefficient (of 2) which is essential to the analogy 
in every possible version; 
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b. the assumption of isentropic gas flows; 

c. the treatment of flows in the equatorial plane without regard for the interaction with the 
general three-dimensional flow field outside of this plane; this shortcoming manifests itself in 
neglecting the term dWdz in the three-dimensional continuity equation; 

d. the omission of radiation pressure, turbulence, viscous friction, heat conduction, diffusion, 
and other real gas effects; 

e. experimental difficulties, e.g., the ones discussed earlier. 

Even though this seems to be a formidable accumulation of shortcomings of the method presented 
in this paper, the approach may be expected to yield more reliable and realistic results about suf- 
ficiently dense atmospheres than the theoretical methods employed by other authors. 


Goddard Space Flight Center 
National Aeronautics and Space Administration 
Greenbelt, Maryland, December 28, 1967 
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Appendix A 

Symbol List 

A characteristic constant defined in Equation 26 
B characteristic constant defined in Equation 32 
c pg specific heat at constant pressure 

F. surface function defined in connection with Equation 10 
3 Froude-number 
g gravitational acceleration 

H w depth of water rotating as a solid body; defined in Equation 13 
H max maximum altitude of atmosphere 
M mass of celestial body 
p pressure 

p 00 atmospheric pressure at the free surface of water 
7 distance from center of rotation 
r 0 reference length 
R g gas constant 
7 time 

T absolute temperature 

u, v, w velocity components parallel to the coordinate axes 
x, y, z Cartesian coordinates 

r universal gravitational constant 

e small dimensionless length, defined in connection with Equation 34 
n w disturbance amplitude of water depth 
k polytropic exponent of a gas 
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wavelength corresponding to 


K 

w 

(jl mass ratio 
p density 

V nondimensional gravitational potential 
w gravitational potential 
co constant angular velocity 

Subscripts 

reference quantities or quantities taken 
at the surface of the earth 

bottom of container of water 

free surface of water 

free surface of water in case of rotation 
of water as a solid body 

All dimensional quantities are marked by a bar above the letters. 


i gas and water 
g gas 
w water 


a, ft first and second celestial 
body, respectively 


1 

2 
II 


llllllllllllllllllllllll 
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